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Abstract. We study rcgularizations of Schwartz distributions on a complete 
Ricmannian manifold M. These approximations are based on families of 
smoothing operators obtained from the solution operator to the wave equation 
on M derived from the metric Laplacian. The resulting global rcgularization 
processes are optimal in the sense that they preserve the microlocal structure 
of distributions, commute with isometries and provide sheaf embeddings into 
algebras of generalized functions on M. 
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1. Introduction 

In this paper we introduce a method of regularizing distributions on a smooth 
manifold by nets of smooth functions such that the approximating nets themselves 
retain a maximal amount of analytic information about the distribution. In par- 
ticular, the analytical properties of interest include the support of the distribution, 
its microlocal singularities (wavefront set), its Sobolev regularity and its behavior 
(pull-back) under certain diffeomorphisms. We shall first abstractly describe what 
properties such an approximation should have and then construct such approxi- 
mations by a suitable choice of smoothing process. The latter is obtained using 
functional calculus for the solution operator of the wave equation for the metric 
Laplacian (we note that the set of analytical properties to be preserved excludes, 
e.g., smoothing via the heat kernel as a possible approximation procedure). 

In order to give a precise formulation of the requirements to be imposed on our 
smoothing processes we need a conceptual framework that allows to assign geo- 
metrical and analytical properties like those mentioned above to regularizations, 
that is, to nets of smooth functions. Such a framework is in fact available in the 
theory of algebras of generalized functions ([3, 4, 27, 25, 15]), which therefore will 
provide the underlying language for our approach. The basic idea in this theory is 
to express analytical properties of distributions as asymptotic estimates in terms of 
a rcgularization parameter e. Up to now, there is a certain dichotomy in the theory 
of algebras of generalized functions. On the one hand, so-called full Colombeau 
algebras allow a canonical embedding of the space of Schwartz distributions on dif- 
ferentiable manifolds ([14, 16]), but their elements do not depend on a single real 
rcgularization parameter s. Instead, such generalized functions are smooth maps on 
certain spaces of test functions and require a rather involved asymptotic. So-called 
special Colombeau algebras, on the other hand, are modelled directly as quotients 
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of certain powers of the space of smooth functions, hence allow for a more straight- 
forward modelling of singularities. A rich geometric and analytic theory is available 
for such algebras (e.g., [27, 25, 5, 7, 13, 19, 23, 20, 12]). The drawback here is that 
there is no canonical embedding of distributions into such algebras. Up to recently, 
only 'non-geometric' embeddings, based e.g. on de Rham regularizations (basically 
through convolution with a mollificr in charts) were available, cf. [7, 15]. In [G], 
however, a new approach to embedding distributions into special Colombeau alge- 
bras was put forward, namely a geometric embedding of distributions on compact 
manifolds without boundary based on functional calculus of the Laplacian. In the 
present paper we follow this general philosophy to produce geometrical embeddings 
for general complete Riemannian manifolds. A main new ingredient here is that we 
employ the solution operator for a certain initial value problem of the wave equa- 
tion for our regularization processes. We obtain a set of optimal properties for such 
embeddings. In particular, they commute with isometrics, respect the functional 
calculus of the Laplacian, and preserve the microlocal structure of distributions. 

The paper is organized as follows: in the remainder of this introduction we fix 
some notations and terminology. Section 2 collects a number of results on wave 
equations on complete Riemannian manifolds. These preparations are then used 
in section 3 to construct optimal regularization processes and use these to obtain 
geometrical embeddings of Schwartz distributions into special Colombeau algebras. 
Finally, section 4 shows how to extend our approach in various directions. On the 
one hand, we demonstrate how to adapt the construction to obtain embeddings 
for distributional sections of vector bundles. On the other hand, we specify the 
main properties of the Laplacian that were used to obtain optimal regularization 
processes in section 3 and show that a wide class of differential operators allows to 
obtain analogous regularization processes. 

Throughout this paper M will denote an oricntablc complete Riemannian man- 
ifold of dimension n with Riemannian metric g. The space T>'(M) of Schwartz 
distributions on M is defined as the dual of the space i7"(M) of compactly sup- 
ported n-forms on M . We write V{M) for the space of smooth compactly supported 
functions on M . Since M is orientable and Riemannian, we may identify T>(M) 
with f2"(Ajf) via /*—>/• dg, with dg the Riemannian volume form induced by g. In 
this sense, D'(M) is in fact the dual space of D(M). We consider L\ oc {M) (hence 
in particular C°°(M)) a subspacc of M via / i— >■ [<p i— > J M ftpdg}. If E is a vector 
bundle over M then T>'(M : E), the space of E- valued distributions on M is given 
by V{M : E) = V (M) ® C °°{M) T°°(M : E), with T°°(M : E) the space of smooth 
sections of E (cf., e.g., [15], 3.1 for details). The wavefront set of a distribution 
w G V(M) is denoted by WF(u>). 

We now turn to notations from the theory of algebras of generalized functions, 
where we basically adopt the terminology from [15, 11]. Given E a locally convex 
(Hausdorff) topological vector space, one can associate to £ a space Ge of general- 
ized functions as follows. Let / be the interval (0, 1]. Define the smooth moderate 
nets in E to be smooth maps (in the sense of [22]) 

I — >• E mu E 

such that for all continuous semi-norms p on E there exists an integer N such that 



(1) 



\p(u e )\ = 0(e N ) as 0. 
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Here as usual by f(e) = 0(g(e)) as e — > we mean there exists an eq > and 
a constant C > such that f(e) < Cg(e) for e < eq. We denote the set of all 
moderate smooth nets in E by Me- Similarly wc can define the negligible nets to 
be the smooth maps u e such that (1) holds for all continuous seminorms p on E 
and all N. Wc shall denote the set of all smooth negligible nets by Me- 

The space of generalized functions based on E is then defined to be the quotient, 

Q E := Me/Me- 

If E is a locally convex algebra then is an algebra as well. One notes that in 
defining M e and Me it suffices to restrict to any family of seminorms that generate 
the locally convex topology on E. If (u e ) is a moderate net in Me then the element 
it represents in the quotient Ge will be written as [(« E )]. 

When E = C°° (M) is the algebra of smooth functions on a manifold M then we 
write M C oo [M) = £ M {M), M C oo {M) = Af(M), and Q(M) := Q C oo {M )- G(M) is the 
standard (special) Colombcau algebra of generalized functions on M ([3, 7, 15]). 
For E — C the space Qc inherits a ring structure from C and we call it the space of 
generalized numbers and denote it by C. Every space Ge is naturally a C- module, 
and hence is often referred to as the C-modulc associated with E ([11]). 

We recall the functoriality of the above construction. If (f> : E —> F is a continuous 
linear map between locally convex spaces E and F then there is a natural induced 
map </>* : Qe —> Qf defined on the representatives as </>*([(w e )]) = [(4>(u e ))]. For 
example any smooth map between two manifolds / : M — > N gives rise to a pullback 
map /* : Q{N) — > Q(M). As a consequence we can define a presheaf of algebras on 
M by assigning to any open set U C M the space G(U). The restriction maps are 
given by the pull back under inclusions, that is if i : U <->■ V is an inclusion of open 
sets then i* : Q{V) — > Q(U) is the restriction map. This presheaf is in fact a fine 
sheaf. Thus in particular we can define the support of a global section u e Q{M) 
as usual to be the complement of the biggest open subset of M on which u restricts 
to 0. In a similar fashion if E — >• M is a (complex) vector bundle then we obtain a 
sheaf of C-modules defined as Q(M : E) := G r ~(M:E) = G{M)® C °°{M) T°°{M : E). 

For any locally convex space E we can also define a subspace <7|? of regular 
elements of Qe ■ These are all elements in Qe such that there exists an integer N so 
that (1) holds independently of the scminorm p chosen. Again wc shall denote by 
Q°°(M) the algebra GeLr M y The algebra Q°°(M) provides the regularity features 
for the analysis of generalized functions in Q(M) in the same way that C°°(M) 
provides these features in V (M) ([27, 18, 5, 13]). For instance: 

(a) Singular support: For u £ Q(M) the singular support is defined as the comple- 
ment of the largest open set U on which the restriction u\u is in Q°°(U). 

(b) Wavefront set: Let fi C 1™ be open. A generalized function u £ G(Ci) is called 
£°°-microlocally regular at (xo,£o) G \ if there exists some <p £ X>(fi) 
with ip(xo) = 1 and a conic neighborhood V C M™\0 of £o such that the Fourier 
transform iF(tpu) is rapidly decreasing in T, i.e., there exists N such that for 
all /, 

(2) sup(l + |e|)'|(^ £ ) A (C)| = 0(e- N ) (e -+ 0). 

The generalized wave front set of u, WF 9 (w), is the complement of the set of 
points (xq,^o) where u is <7°°-microlocally regular. 
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An alternative description of WF g (u) is as follows ([13]): Let P be an order 
classical pseudodiffcrential operator and let char(P) C T*M be the charac- 
teristic set of P, that is the 0-set in T*M \ of its principal symbol. Then for 
u G Q{M), 

WF 9 (u) = P| char(P) Pe^(M). 

(c) Hypoellipticity: An operator P is said to be C/°°-hypoclliptic if for every U C M 
open and every u G G(U), 

Pueg^iUj^ueg^iU). 

General references for microlocal analysis in algebras of generalized functions are 

[25, 5, 18, 19, 12, 13]. 

2. The wave equation on a complete Riemannian manifold 

In our approach, optimal regularization processes on complete Riemannian man- 
ifolds will be based on the solution operator for the wave equation. To allow for a 
smooth presentation, in the present section we therefore collect some basic proper- 
ties of solutions of the wave equation in this global setting. 

Let (M, g) be an oriented, connected complete Riemannian manifold (without 
boundary) of dimension n and denote by A the Laplace operator on M. The 
Riemannian metric g induces a volume form dg on M, and we will denote the 
corresponding L 2 -norm by || ||l 2 (a/)- On differential forms, the corresponding inner 
product is given by (a, 0) := J a * (3 = J a A *(3. 

Let d be the exterior differential on the space W(M) of differential forms on 
M and denote by * the Hodge star operator. Then the codifferential 5 on f2*(M) 
is defined, for any fc-form a, by 5a = (— \y ik + n + 1 * d * a. Finally, the Laplace 
operator on Q* (M) is defined by A := (d + S) 2 = do S + S o d. This sign convention 
renders A a positive operator on L 2 (M) (cf. [10]), and for any smooth function u 
in particular Au = — divgradu. 

The operators d, S and A are unbounded on the Hilbert space L 2 (M : A*M). 
The natural domain of d is given by Dom(d) := {a <E fi*(M) | ||a|| , ||da|| < oo}, 
and analogously for 5. This fixes the natural domain of A to be 

Dom(A) := {a G Dom(d) | da G Dom(5)} n{a£ Dom(<5) | 5a G Dom(d)}. 

We will mainly be interested in the restriction of A to L 2 (M,dg), which is an 
unbounded essentially self-adjoint operator with dense domain (cf. [10]). 

We consider the following initial value problem for the wave equation on M (or, 
strictly speaking, on R x M): 

(3) ( ^+A)u = 

(4) u(0,x)=«o(s) ^-u{0,x) = 

os 

Since g is complete, the Laplace operator is self-adjoint and the above wave equation 
has a unique (mild, hence distributional) solution in C(R, L 2 (M)) for all uo in 
L 2 (M). By functional calculus this solution can be written as cos(sy/A)uQ. 
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Remark 2.1. We briefly sketch a proof to the existence and uniqueness result for 
the above wave equation. Since A is a positive self-adjoint operator, we may equiv- 
alcntly rewrite (3)-(4) as a first order initial value problem on the Hilbert space 
H :=L 2 {M)®L 2 (M): 



d I u\ ( — n/A / 



ds \v J \ ivA J\y 

u(0, . ) = M v(0, . )vq := iV~Ku 



We set 



/0 A , f-y/K 
A) := n n := „ ^ w := 



Uq 

0/ — "\ o VA/ ~ l "~Wo 



By the theory of unitary semi groups (e.g., [28]), %A\ generates a strongly continuous 
unitary group U(s) = exp(isA\). Since Aq is bounded, A := A^ + iAi also generates 
a strongly continuous semigroup. Consequently, the above initial value problem 
with wq € Dom(A) is uniquely solvable. More explicitly, for wq in the dense 
subspacc D°° := H^Lo Dom(A fc ), the power series expansion of exp(s^4)u;o readily 
shows that u(s) = cos(sv / A)tto on a dense subspace, hence in fact for all uo € L 2 (M) 
(cf. also [2, 30]). 



For a given even Schwartz function F € J^(R) the operator F(vA) can be 
defined by functional calculus for essentially self adjoint unbounded operators. We 
show that we have the following alternative description of F(\/~A) by inverse Fourier 
transform: 

(5) F(VA) = — / F(s) 003(5^) ds. 

In fact, by standard estimates in the functional calculus this identity holds pointwise 
on L 2 (M). Moreover, since 



||F( S )cos(sVA)|| <|F(«)|||cos(«.)IU- <\F{s)\ 

and F £ L 1 (R), the integral in (5) exists as a Bochner integral in B(L 2 (M)). 

The description of functional calculus using (5) allows to estimate the support 
of the kernel of the operator F(vA) based on the finite speed of propagation for 
the operator cos(sv / A). To see this, we first describe an alternative approach to 
obtaining the solution operator to (3)— (4). Consider the first order differential 
operator D := d + S on the space f2*(M) of differential forms on M. The symbol 
an of D is given by ao(x,^) — £ A . —t£, where denotes interior differentiation 
along the vector field metrically equivalent to the one-form £ (cf., e.g., [26], 10.1.22). 
Therefore, 

<r D (x,0 2 = (tA . -* e ) 2 = -iien 2 id 

(again by [26], 10.1.22). We conclude that the speed of propagation of D, defined by 
c D := sup{||cr.D(a;,£)|| | x € M, £ e T X M*, j|£|| = 1} is c D = 1. Since D is symmetric 
and of finite propagation speed, it is essentially self-adjoint ([17], 10.2.11). 

Based on these facts, an explicit bound on the speed of propagation for the 
support of e lsD u is given by the following result (see [17], 10.5.4): 

Proposition 2.2. Let u € L 2 (M : A*M) and denote by d g the distance function 
induced by g. Then supp(e lsD u) C B\ s \(supp(u)) := {x € M | d g (x, supp{u)) < \s\}. 
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Thus the same property holds for the bounded operator 

cos(sD) = -(e isD + e~ isD ) 

By functional calculus, for any u <E L 2 (M), cos(sD)u solves the initial value problem 
(3)-(4). By uniqueness, therefore, cos(sD) = cos(sy/~A) on L 2 (M) and the above 
considerations apply to our solution operator. 

Next we provide some estimates that will repeatedly be useful in our further 
study of operators of the form F(y/~A). 

Lemma 2.3. Let F G 5(K) be even. Then for any compactly supported smooth 
function u, 

1 f°° 

\\F(VA)u\\ L 2 {M) < ||u|| L 2 (M) - / \F(s)\ds. 

f Jo 

Moreover, for any positive integers k, I, 

\\A k F(VA)A l u\\ LHM) < \\u\\ L2{M) - / \F^ k+2l \s)\ds. 

71 Jo 



Proof. The first estimate follows from (5): Since the operator cos(sv A) is unitary 
we have 

\\F(VK)u\\ = II— / F(s)cos(sVA)uds\\ 

<||u||i r\F(s)\ds. 
K Jo 

Concerning the second inequality, note that by functional calculus we have 

A k F(VA)A l u = (t 2k+21 F{t))(y/A)u 
1 1"°° 

= _ / (t 2k+2l F(t)) A cos(sVA)uds 
2tt J_oo 



1 f°° 

— / F<- 2k+2l \s)cos(sVA)uds, 
2tt J^oo 



2tt 

from which the claim follows as above. □ 
In particular, the above proof shows that 

FW(VA) = A k F(VA). 
For any s £ M. and any u £ V{M) we set 

H s :=||(l + A) s / 2 U || i2(M) 

The Sobolev space of order s is the completion of T>{M) with respect to this norm. 
We set H°°(M) := f] seR H s (M) and denote by H* p (M) the space of compactly 
supported elements of H S (M). 

The following result will be essential for our approach to regularizing distribu- 
tions on complete Riemannian manifolds. For the notion of (properly supported) 
smoothing (or regularizing) operator we refer to [I], ch. 1.4. 

Proposition 2.4. Let F e S(R) be even. Then 

(i) The operator F(\/~A) : T)'(M) — > C°°(M) is a smoothing operator. 
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(ii) Let c > and let <f) c G T>(R) be such that supp((f> c ) [—2c, 2c]. Then 

1 f°° 

T(a/A) := — / (£ c (s)F(s) cos(s-\/A) ds, 
2tt J-oo 

zs a properly supported smoothing operator. 

Proof, (i) This is a consequence of ellipticity of the Laplace operator A. For brevity, 
we set A := F(y/A). Given any ip e V(M x M), let 

(A» := / (A<p(x, .))(x)dx. 
Jm 

Then for all fa, fa e V(M), 

(K,fa®fa) = / fa{x)Afa{x)dx = (A fa, fa), 

so K is the distributional kernel of A : V(M) -> C°°(M) C V'{M). We have to 
show that A" e C°°(M x M). By Lemma 2.3, for all k e N a and all ip € 2?(M), 
II^IU < C k \\i})\\ L 2 {M) , hence 

KJf,^iO^>l< / \fa{x)\\Afa{x)\dx< C4fa\\ L 2 (M) \\fa\\ L 2 [M] 
Jm 

Since V(M) (g) P(M) is dense in A 2 (M x A/) this implies that if G L 2 (M x M). 
Now set A fe! ; := A l AA k . Then the kernel of A ky i is given by K k j := AyA l x K, since 

{K k>l ,fa®fa) = (K,A l fa®A k fa) = (AA k fa,A l fa) = {A k jfa,fa). 

As above it follows that K k: i € L 2 (M x M) for all fc, Z, hence by elliptic regularity 
for A x ® 1 + 1 Aj, it follows that AT is smooth. 

(ii) T{yf~A) is a smoothing operator by (i). To establish proper support, by [9], 
Prop. 8.12 we have to show: 

(a) VA' CC M 3L CC M such that u e V(K) T(y/A)u £ V(L). 

(b) VA' CC M 3L CC M such that u = on L => T(\fK)u = on K. 

Both (a) and (b) follow from the finite speed of propagation of cos(y / A) which 
implies that there exists some C > such that for any u E L 2 (M), the support of 
T(y/~K)u is contained in a ball of radius C around supp(u) (Prop. 2.2). The result 
therefore follows from the properness of the complete metric g. □ 

3. Embeddings 

In this section we will employ the smoothing operators developed in Section 2 to 
construct optimal embeddings of the space V(M) of distributions on a complete 
Riemannian manifold M into the algebra Q(M) of generalized functions on M. 

A set X C M x M is called proper if the restriction of the projections on both 
factors TTj : X M : j = 1,2 are proper maps. Let &~™ p (M) be the space 
of all operators T : C°°(M) — > T>(M) with smooth kernels with proper support 
in M x M . By a regularization process we mean a net T e of properly supported 
smoothing operators which provides an approximate identity on compactly sup- 
ported distributions. More precisely, we shall be interested in rapidly converging 
regularization processes of the following kind: 

Definition 3.1. A parametrized family (T £ ) e ei of properly supported smoothing 
operators is called an optimal regularization process if 
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(A) The regularization of any compactly supported distribution is of moderate 
growth. That is, for any continuous semi-norm p on C°°(M) and any distribu- 
tion w G S'(M) there exists an integer N such that 

P {T e w) = 0(e N ) (e^O), 

i.e., (T £ w) G £ M {M). 

(B) The net (T e ) is an approximate identity: for each compactly supported distri- 
bution w G £"(M) 

\imT £ w = w mV(M). 

(C) Preservation of supports: For any w G £"(M), supp(i/;) equals the Q(M)- 
support supp[(T £ w)] of the class [(T £ w)}. 

(D) If u G T>(M) is a smooth compactly supported function on M then for all 
continuous semi- norms p on C°°(M) and given any integer m, 

p(T E u-u) = 0(e m ), 

i.e., (T e u-u) eAf(M). 

(E) Preservation of wavefront sets: Setting lt '■ w i— )■ [(T e w)], for any u> G <§'{M) 
we have 

WF(iu) =WF 9 (i T (w)). 
Given an optimal regularization process (T E ), we obtain a linear embedding 

l t : S'{M) -> g(M) 
= [(r e w)] 

(by (A) and (B)). By (C), t T extends to an embedding of V(M) into g(M) which 
preserves supports. More precisely, there is a unique sheaf morphism on T>'(M) (also 
denoted by lt) which extends lt : S"(M) — > G(M). lt is a linear embedding that 
commutes with restrictions. Details on how to extend lt from S"(M) to T>'(M) can 
be found in [7], Sec. 2 or in [15], 1.2 (although carried out for special cases of optimal 
embeddings in these references, the arguments given there, entirely sheaf-theoretic 
in nature, carry over to the general situation studied here). (D) implies that lt 
renders C°°(M) a faithful subalgcbra of Q(M). Finally, (E) secures preservation 
of wavefront sets and, therefore, of singular supports under this embedding. In 
particular, precisely the distributions that map into the subalgcbra Q°°(M) under 
lt are smooth: 

l t (V(M)) n g°°{M) = l t {c°°{m)). 

First, we provide some examples of optimal regularization processes: 

Example 3.2. Let A be the Laplace operator associated to a closed Riemannian 
manifold M. Let F G J^(R) be a Schwartz function on the reals such that F 
is identically 1 near the origin. Let F e (x) :— F(ex). Then by applying standard 
functional calculus, F £ (A) is an optimal regularizing process: For a closed manifold, 
Wcyl's estimates on the spectrum of the Laplacian provide asymptotic bounds for 
the spectral counting function 

^a(A) = #{Afc| Afc < A}. 

In fact, for m = dim(M) we have 

A A (A) ^ Af. 

M j (47r)tr(m/2 + l) 
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Essentially, this suffices to obtain all the estimates in Definition 3.1 (preservation of 
wavefront sets follows as in Th. 3.10 below). In addition, F £ (A) is invariant under 
isometries. We refer to [G] for details. 

Example 3.3. As in the original construction of Colombeau (cf., e.g., [3, 4, 15]) 
an optimal regularization process can be constructed from a mollifier p £ S^(W l ) 
satisfying the following conditions: 

(6) / p{x)dx = l { x a p(x)dx = aeN'j. 

Then the net of functions p £ (x) :— 4rP(f ) i s a delta net. Convolution with such a 
delta net provides an example of an optimal regularization process. For example, 
estimate (D) from Definition 3.1 can be established in this setting using Taylor's 
theorem and the moment conditions (6) imposed on p. Concerning (E) from Dcf. 
3.1, see [25, 18]. An important characteristic of these approximate units is their 
equi variance with respect to the Euclidean translations. 

For (M, g) a complete Riemannian manifold with Laplacian A, F £ S(K.) even, 
and </> c as in Prop. 2.4, we additionally suppose that F equals 1 in a neighborhood 
of and that C is even. For any e £ I we set F £ (s) :— F(es), and 

1 f°° 

(7) T £ (VA):=— ^ c ( s )(F £ ) A ( S )cos( S VA)ds. 

27T J-oo 

We will show that the family of smoothing operators (T e (VA)) ee / is an optimal 
regularization process in the sense of Def. 3.1. We note that each T E (VA) is a 
properly supported smoothing operator by Prop. 2.4 (ii). Turning first to Def. 3.1 
(A), we have: 

Proposition 3.4. Let u £ S"(M). Then (T e (y/A)u) G £ M {M). 

Proof. It follows immediately from (7) that (e,x) i-> (T e (y/~A)w)(x) is smooth, so 
it remains to establish the moderateness estimates for the net (T s (y/A)w). Since 
w € S"(M) there exists some s £l with w £ H°°{M). By Prop. 2.2, there exists 
some fixed compact set K CC M such that supp(T e (v / A)u') C K for all e £ I. 
Moreover, Prop. 2.4 implies that each T e {\f~K)w is in H™(M). By the local Sobolev 
embedding theorem it therefore suffices to show that for each s £ R there exists 
some N £ N such that 

\\T £ (^A)w\\ s = 0(e- N ). 
In fact (by enlarging s if necessary) we may assume in addition that / := £ N. 
Let u be the unique element of L 2 (M) such that w = (1 + A)~ s °/ 2 u. Then 

\\T e (VA)w\\. = ||(l + A) s / 2 T e (VA)(l + A)- s «/ 2 w )|| L2(M) 

= ||(l + A)'T e (>/A)ii|| La(a0 

Now write 4> c = (V'c) A for some tp c £ £>(]R). Then 

T £ (VA)u= / (ip c * F £ ) A (t)cos(tVA)udt. 

J -2c 
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Since ip c * F s is even, Lemma 2.3 implies that 

||A^T e (v / A>|U2 (M) < || U || L2(M) - / \[(^ c *F e ) A ]^(t)\dt 

71 Jo 

From this, we finally obtain 

(8) \\AiT £ (VA)u\\ LHM) < \\u\\ LHM) ± I* \Mc{t)\p(i))\&Xt)\dt 



7T 

= Oie- 21 - 1 ) 

for o < j < i. a 

We may use the method of proof of Prop. 3.4 to show that the embedding it is 
in fact independent of the particular choice of the cut-off function (j> c : 

Lemma 3.5. Suppose that 0^ and (f> 2 2 are cut-off functions as in (7) and denote 
the corresponding regularization operators by T* (\/A) and T 2 (\/~A), respectively. 
Then for any w £ £"(M), l t i(w) = l T 2(w). 

Proof. Set 4> <Pci ~ ^C2- Using the assumptions and notations from the proof of 
Prop. 3.4, we have to estimate 

\\AJ(T^VA)u-T?(VA)u)\\ LHM) . 

Since 

T^(VA)u-T*(VE)u = — / <Xs)(F 6 ) A ( s )cos(sv / A)uds, 
analogous to (8) the L 2 -norm of this expression is bounded by 
\Hlhm)- / \[m)-F(~))}M(t)\dt. 



7T _ 

A typical term to estimate therefore is 

\ij><p)(er)Fto(r)\dr ( P ,qeN ) 

Since F = 1 near 0, all higher moments /_ r k F(r) dr of F (k > 1) vanish. More- 
over, <fi is identically zero in a neighborhood of 0, so Taylor expansion of 4>( p >{er) 
around zero implies that the above integral is of order e m for any m £ N. □ 

Next, we establish suitable convergence of lt(u>) to w. 

Proposition 3.6. Let w £ S'{M). Then T E (s/A)w -> w in V'{M). 

Proof. We may write w = (1 + A) k u for some u £ L 2 (M) and some k £ No. By 



Prop. 2.3 it therefore suffices to show that T e (\f~R)u — u in L 2 (A1) in order to 
ensure that T e (y/A)w -> w in H k (M) and hence in V(M). Now 



\\T E (vA)u - u\\ L 2 {M) 

/OO 
\\(f> c (er)F(r)(cos(erV~A)u - u)\\ L *( M ) dr -> 

by dominated convergence. □ 

This settles Def. 3.1, (B). As was remarked after Def. 3.1, we thereby obtain a 
linear embedding lt of S"(M) into Q(M). Our next result establishes preservation 
of supports under lt- 
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Proposition 3.7. For any w £ S"(M), supp(w) = supp(lt{w)). 

Proof. Let x £ M \ supp(ui) and choose a compact neighborhood K of x such that 
K n supp(ui) = 0. Suppose first that w is continuous. Then 



Lr{w) e (x) = I 4> c (sr)F(r) cosferV A)w(x) dr. 

J —oo 

We split this integral in one part over |r| < 2c/y/e and a second part where 
\r\ > 2c/y/e. For the first part we note that by Prop. 2.2, supp(cos(er-\/A)u>) C 
5 e | r |(supp(u;)) C B^2c( SVL PP( w )) f° r au M < 2c/ v^- Hence for small e, the sup- 
port of the first term lies in the complement of K . 

To estimate the second integral we proceed as in the proof of Prop. 3.4: observe 
first that by Prop. 2.2, the support of cos(£T-\/A) is contained in a single compact 
set for all e and all r in the domain of integration. It therefore suffices to estimate, 
for each j £ Nq- 



</> c (£r)-F(r)(cos(erv A)A J w) dr\\ L 2 {M) 

M>2c/y? 



<\\w\\ 2j \F(r)\dr = 0(e m ) 

J\r\>2c/^ 

for each m. Summing up it follows that x does not lie in the support of lt(w) 
in Q(M). In the general case where w is not necessarily continuous we can write 
w = (1 + A) k v for some continuous v and some k £ Nq, so the above argument 
readily carries over. 

Conversely, let x £ supp(w) and suppose that there exists a neighborhood U of 
x such that Lt{w)\u = in Q{M). Pick some cp £ T>(U) such that (w, <p) ^ 0. Then 
\{lt(w) e , <f)\ = 0(e m ) for each m but (lt(w) —w,<p) — > by Prop. 3.6, so we arrive 
at a contradiction. □ 



Proposition 3.8. Let u £ V(M). Then (T e (VA> - u) eeI £ Af(M). 

Proof. By the local Sobolev embedding theorem, it suffices to show that for all 
J G No 

a(e,j) := ||A''(T e (\/A)u-tO||i> (Jl0 - 0(e m ) 
for each m £ N. Due to our assumptions on F and (f> c , a(e,j) equals 



/ (F £ ) A (i)(</> c (t)cos(tv / A) w -A^)|| L3(M) 



F(r)4> c (er)(cos(erV A)A 3 u - A J u) dr\\ L 2 (Af) 



By Taylor expansion, for any m £ N there exists some C m such that 

z l r l 



m— 1 

(er) cos(er\/A)A J u = A 3 u + ^-^A^^u + i? m (r, e)A 1+m l 2 u 

i=i 



where aj £ M. and R m is globally bounded by C m e r ' 
Since all higher moments of F vanish, 



a(e,j) < || / F(r)R m (r,e)Ai+ m / 2 udr\\ L 2 {M) 

J —GO 

< C m ||^|| il(R) ||A^™/ 2 u|| i2(M)£ m , 
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as claimed. □ 

The following important invariancc properties of the embedding ij- follow imme- 
diately from (7). 

Proposition 3.9. 

(i) Let f : M — > M be an isometry. Then for any u G 'D'(M), iT{f*u) = 

f*L T ( U ). 

(ii) If is a pseudodifferential operator commuting with A, then ^ commutes 
with it- 

Turning now to the singularity structure of distributions and their embedded 
images, we will show that the embedding lt preserves the wavefront set of distri- 
butions, i.e., that (E) from Def. 3.1 is satisfied. 

Theorem 3.10. Let w G V'(M). Then WF(w) = WF g (i T {w)). 

Proof. We first note that the notion of wavefront set (both distributional and gen- 
eralized) is local in nature. Moreover, by finite propagation speed of the solution 
operator cos(sy/~K) (Prop. 2.2), we may choose the cutoff function <f> c in such a way 
that given some local chart (ip, U) and any w G S"(U), each T e (\fK)w is supported 
in U as well (this particular choice of <fi c does not affect Lt by Lemma 3.5). By 
unique solvability of the wave equation (3)-(4) we may therefore use (7) on tp(U) 
(with metric ip*g), thereby effectively transferring the problem to R'\ 

Suppose first that (x ,£o) G (T*M \ {0}) \ WF s (w). By (2) this means there 
exists some conic neighborhood T of £o in T*M \ {0}, some N G No and some 
ip G V{R n ) with ip(x ) = 1 such that for all I G N and all £ G T, 



(9) \(ip(T s (VA)w)) A (0\(l + ICI) 1 = 0(e 

For a suitable k G No we may write w = (1 + A) k u, with u G H 2 (R n ). Since A 
is elliptic, WF(w) = WF(ti) and WF 9 (i T (w)) = WF 9 (i T (ii)) (by [13], Th. 4.1 and 
Prop. 3.9 (ii)), so we may without loss of generality assume that w G H 2 (W l ). We 
have to estimate 



K^) A (oi < m(VA) W ~w)^(o\ + \(MVA)wr(o\ 

(10) < \\<p(T e (VA)w - w)\\ L i + |(^T e (v / AH A (0|. 

Now 



ip(x)(T e (V A)w(x) — w(x))\ dx 
— J J Wi^W&c^r) cos(er\ r A)w(x) — w(x)\ dx\F{r)\ dr 
<C(ip) / \\<l> c (er)cos(er\ r A)w -w\\ L 2\F(r)\dr 



We use functional calculus to bound this term (cf., e.g., [2!)], Th. VIII. 4). Let 
U : L 2 (W l ,g) L 2 (fl,dfj,) be a unitary isomorphism transforming A into the 
multiplication operator h n- fh (for some fixed / G L 2 (il,dfi)). Then setting 
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a(e, r) := \\4> c (er) cos(er\/A)w — w\\ L 2 we find 

a(e,r) 2 = J \^ c {er) C os{er^J(u)) - l\ 2 \{Uw){u)\ 2 d^Lo) 

< Ce 2 r 2 J f{uf\{Uw){u)\ 2 dft{w) 

Since w G Dom(A), / • (Uw) G L 2 (fl,dfi), hence the integral in this last expression 
is finite. It follows that a(e,r) < Ce\r\. From (10) and these calculations we 
therefore conclude that for some C — C(ip, F), 

(11) |(^) A (0|<C £ + |^T e (VAH A (e)|. 

We now show that for any m G No, |£| N+1 |(<pw) A (C)l i s bounded on T, thereby 
demonstrating that (£o,£o) ^ WF(w). 

Suppose to the contrary that there exists some m G No and a sequence £j G 
T with -> oo such that l^jl^r |(<£tu) A (£j)| -> oo as j ^ oo. Then £j := 
l^ r™^ 1 — * 0) an d using (11), we obtain 

\z j \^ k \(<pu) A &)\ = #ie,f m iw A fe)i 

< C £ f +1 |e j | 2m +ef |0| 2m |(^(^) W ) A (OI 

By (9), however, the right hand side of this inequality is globally bounded, a con- 
tradiction. 

Conversely, suppose that (a;o,£o) ^ WF(w). We have to show that (xo,£o) 
WF g (tr(u>)). As above, we may without loss of generality suppose that w G L 2 (M). 
Pick some open neighborhood U of xq and some conic neighborhood Ti of £o hi 
E™ \ such that (U x Ti) n WF(w) = 0. Let us suppose for the moment that we 
already know that, setting u(s, x) := cos(s\/A)w, we have 

(12) 3 So >0:(C/xr 1 )n{(x,O|3 S , | s| < s 3r : (s,x;t, £) G WF(ii)} = 0. 

Then, given x S 2?(t/) and ^ G T>((— s , s )), for each Z G N there exists some 
Ci > such that 

|((«/®x)tO A (r,OI <C,(l + |r| + |e|r' (rGR, £ G i\). 
Thus for / > n and |s| < so we obtain 

K«)(x-«(«,.)) A (0I = I^ s ((MxW a K-t))I 

In addition, we now choose c such that 2c < so (which is possible by Lemma 3.5) 
and v G T>((—sq, sq)) such that v = 1 on supp0 c . Then for £ G Ii, 

|(x-^He) A (OI - I / <f>c(s)(F £ ) A (s)is(s) I e-* x x {x)u{s,x)dxds\ 

= oca + 

Thus, (xo,£o) ^ WF g (tT(w)), as claimed. 

It remains to establish (12). To this end, denote by (3 the bicharacteristic flow 
on T*(R x M) corresponding to d 2 + A. Since u is the solution to (3)-(4) with 
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u = w, by [8], p 118, WF(ii) C C o WF(w), where 

Co = {((s,x;t,£), (x Ao)) I 3r, t e K : (s,x',t,£) = 0(r, (0,x ,t q ,^ o )) 

A -t 2 +s xo (£ ,£o) = 0}. 

/3 is the flow of the Hamiltonian vector field of the symbol — r 2 + g x (£, £), so the 
corresponding system of ODEs reads 

s(r) = — 2r(r) 

x(r) = 2g x(r) {^{r), .) 

f(r) = 

i(r) = -Dg(x(r)Mr),ttr)) 
Denoting by Pi the z-th component of j3, it follows that 

C = {((s,x;t,£), (xo,£o)) I r = r = ± vV^o, £o), s = -2rr ,3r 6 K : 

= G8 a ,&)(r, (0, x ,t , &>))}, 

and, since WF(u) C C o WF(ra), 

WF(«) C {(s,x;r ,e) | r = ±vW&,£o), a = -2rr ,3r e H : 

3(2, D G WF(«)) 3r 6R:(i,{)= (/3 2 ,/3 4 )(r, (0, x, r , f))}. 
By continuity of f3 and the fact that (U x r x ) n WF(ra) = 0, (12) follows. □ 
Summing up, we obtain 

Theorem 3.11. The family (T £ (v / A)) e e/ defined by (7) is an optimal regulariza- 
tion process. The corresponding embedding 

l t : V'(M) -> Q(M) 

l t (u) = [(T e (VA)u)} 

is an infective sheaf morphism that renders C°°(M) a subalgebra of Q(M). lt 
commutes with isometries and pseudodifferential operators that commute with A. 
Moreover, it preserves the singularity structure (wavefront set) of distributions. 

Remark 3.12. We note that while Th. 3.11 is formulated using the language of 
algebras of generalized functions, it can also be used independently of this theory. 
For example, on the level of regularizing nets, preservation of wavefront sets under 
tr means that the wavefront set of a distribution w G V(M) can be read off from 
the asymptotic properties of its regularization (T E (v / A)(w)) via (2), and similar for 
the other properties. 

4. Distributional sections of a vector bundle 

In this section we consider the problem of regularizing distributional sections of a 
vector bundle over a manifold. We shall provide a notion of optimal regularization 
and show that given a differential operator D satisfying two simple conditions, we 
can always obtain such regularizations. 

Let \A\M denote the density bundle over M. Then for E* the dual vector 
bundle of some vector bundle E, the space of distributional sections of E is given 
by V'(M : E) := Yf(M : E* ® |A|A/)'. In particular we have a natural inclusion 
r°°(M : E) -> T>'(M : E). By choosing a trivialization of the density bundle \A\M, 
for example by choosing a Riemannian metric, and by choosing a Hcrmitian inner 
product on E to identify with E* we can (non-canonically) identify T>'(M : E) 
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with r^°(M : E)' . In the sequel we shall assume that we are given a Riemannian 
metric on M and a Hermitian inner product on E. We similarly define the space 
£'{M : E) of compactly supported E- valued distributions. 

By a smoothing operator on E we shall mean an operator defined by a kernel in 
T°°(M x M : End(.E) <g> A R )\ Then if T is a smoothing operator then T : S'(M : 
£) ->■ T°°(M : E). 

Having fixed our notations we shall define optimal rcgularization processes for 
distributional sections analogous to Def. 3.1. 

Definition 4.1. A parametrized family (T e ) e6 j of properly supported smoothing 
operators is called an optimal regularization process if 

(A) The regularization of any compactly supported distributional section s £ 
£'(M : E) is of moderate growth: For any continuous seminorm p on L°°(M : 
E), there exists some integer N such that 

P (T e s) = 0(e N ) (e^O), 

(B) The net (T £ ) is an approximate identity: for each s £ &'{M : E), 

limT £ s = s in V (M : E). 

(C) If u £ r^°(M : E) is a smooth compactly supported section of E then for all 
continuous seminorms p and given any integer m, 

p(T e u-u) = 0{e m ). 

(D) The induced map it ■ £'(M : E) — > Q(M : E) preserves support, singular 
support and the wavefront set. In particular, 

Lt{V'(M : E))n g°°(M : E) = r°°(7\7 : E). 

In the following section we shall describe the precise requirements on a differential 
operator D that would provide us with the functional calculus necessary for the 
construction of an optimal rcgularization. 

4.1. Admissible operators. Let E — > M be a vector bundle over a complete 
Riemannian manifold M provided with a Hermitian inner-product ( }e- We shall 
denote by L 2 (M : E) the completion of the compactly supported sections r^°(M : 
E) with respect to the norm 

\\s\\ := [ (s(x), S (x)) E dx s£T™(M:E). 
Jm 

Definition 4.2. Let D be a symmetric first order differential operator on E. We 
shall assume that 

(1) The operator D has finite speed of propagation, that is the norm of the 
principal symbol over the unit sphere is bounded by a constant Cd, 

C D = Bup{||<7z>(a:,0ll \ x£M, = 1} < oo. 

(2) The operator D is elliptic. 

Such a differential operator D shall be called admissible operator. 

^To be precise, let 'kl^r '■ M X M — > M be the left and right projections on M. Then 
End(E) := n* R (E)* ® tt* l {E) and A R = 7r|j|A|M. 
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As a consequence of the finite speed of propagation, D is essentially self-adjoint. 
Therefore the equation 

d 

(13) —u~iDu u(.,0) = uq, 

has a unique solution for all times t for any initial datum uq G (M : E) . Unique- 
ness follows from energy estimates, while existence is seen by applying functional 
calculus to note that e ztD uo is a solution. 

Furthermore for a Schwartz function F G ,y($L) the Fourier inversion formula 
gives that 

(14) F(D) = — / F(s)e isD ds. 

As already noted, if u is supported in a set L then e ztD u is supported in the ball 
Bc D -t{L). This has the following consequence: 

Lemma 4.3. // F G J5^(R) is a Schwartz function such that F is supported in an 
interval (— c, c) then, for any u G L 2 (M : E), 

supp(F(D)u) C B CD -c(supp(u)). 

On the other hand the ellipticity of D insures that the operator defined by 
applying a Schwartz function F to D is necessarily a smoothing operator. 

Let as before F be an even Schwartz function in y(M.) such that F = 1 near the 
origin. Let F e (x) := F{ex). Our main result in this section is 

Theorem 4.4. Given an admissible differential operator D and a Schwartz function 
F as above the family of operators (F e (D)) £ £j provides an optimal regularization 
process in sense of Def. J^.l. 

The next two subsections provide the arguments for the proof. 

4.2. Weyl's law and functional calculus. In this subsection we shall assume 
that M is compact. In this case any symmetric operator D is essentially self- 
adjoint. In addition the operator D 2 is a positive elliptic operator by assumption 
and hence Weyl's asymptotic formula for eigenvalues gives 

(15) N D2 (X) := #{A, G sp(D 2 )\ A, < A} ~ CA^^. 
Then the following can be obtained by applying (15). 

Lemma 4.5. Let D be an elliptic self-adjoint differential operator of order 1 and 
let M be compact. Then for a Schwartz function F onM. with F = 1 near the origin 
we have: 

(A) Given a smooth section u G Y(M : E) 

\\F e (D)u - u\\ L 2 (M:E) = 0(e m ) for all m G Z. 

(B) If s G" H k {M : E) for every k > t then given any S > 0, \\F 6 (D)s\\ L 2^ M . E ) is 
not 0(e ™ M +*+< 5 ). j n particular, 

l Fe{d) {V'{M : E)) n G°°(M : E) = T°°(M : E)). 

(C) For every distributional section s the regularization (F £ (D)s) is moderate. 
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The proof of the above Lemma can be found in [l i] . This result is precisely due 
to the fact that D 2 is a positive elliptic operator. We still need to prove that the 
microlocal properties hold true for our embeddings F e (D). These turn out to be 
precisely due to the finite speed of propagation of D. 

4.3. Finite speed of propagation and localization. We now return to the gen- 
eral situation where M is a complete Riemannian manifold not necessarily compact. 

Recall that if X is a compact manifold with boundary then one can obtain a 
double of X, denoted here by DX by gluing two copies of X along the boundary 
dX (e.g., [21], VI 5.1). Now if X is a compact manifold with boundary embedded 
in a Riemannian manifold M of the same dimension and if U is an open subset of 
M such that U C interior(X), then one can choose a Riemannian metric on DX so 
that the inclusion j : U DX is an isometry. Furthermore it is clear that given 
any vector bundle E — > M there exists a vector bundle Ex —> DX such that Ex 
restricted to U is canonically isomorphic to E. At the same time there exists a 
symmetric elliptic operator Dx on Ex that matches up with D on U. 

We fix a compactly supported distributional section u € £' (M : E) and a con- 
stant c > 0. Since M is complete the open ball U := Bi c .c D (supp(zt)) is relatively 
compact and is contained in a compact manifold with boundary X C M. Now s 
can be identified with a distributional section of a vector bundle Ex — > DX. 

Proposition 4.6. With assumptions on u, c and F as above, let us further assume 
that the Fourier transform F(s) is supported in an interval (— c, c). Then F(D)u 
and F(Dx)u are both supported in U and 



Proof. Since the operators D and Dx restricted to the open set U coincide, the 
uniqueness of solutions to the equation (13) implies that e zsD u and e lsDx u agree 
for s < Cfj- The statement therefore follows from the Fourier Inversion Formula 



We can now finish the proof of our main result. 

Proof of Theorem 4-4- First we note that given any cutoff function <fi(s) supported 
in an interval (— c, c) such that <f> = 1 near the origin, and any compactly supported 
distributional section u, 



With this observation it is clear that: 

(1) All estimates for Definition 4.1 follow from Lemma 4.5. 

(2) The support of u coincides with the generalized support of [F e (D)u]. This 
implies that the embedding extends to a sheaf morphism ip e :T>' [M : E) — > 



Q{M : E). 

(3) Since wave-front sets are defined locally, our embedding Lp e preserves wave- 
front sets by Th. 3.10. 



F(D)u 



F{D x )u. 



(14). 



□ 




in Q{M : E). 



□ 
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Remark 4.7. From the proof one notices that a second order positive elliptic dif- 
ferential operator A on sections of E also provides us with an optimal embed- 
ding F £ (A) provided that the solution operator to the wave equation (3), namely 
cos(svA) propagates support at a finite speed. Thus in particular if TJ(M) denotes 
the tensor bundle on M and g a complete Riemannian metric on M the induced 
Laplace operator provides an example of such an operator. 

4.4. Isomorphisms between vector bundles. Let : E\ — > E 2 be an isomor- 
phism of Hcrmitian vector bundles (preserving the inner product). Given any ad- 
missible differential operator D\ on sections of E\ the push-forward D 2 := fyD^r 1 is 
also an iso-spectral admissible differential operator. In particular, for any Schwartz 
function F we have 

F(D 2 ) = 4>F{Dx)4>- 1 . 

The extension of 0* : F°°(Af : Ei) F°°(M : E 2 ) to the generalized sections, 
4>* : Q(M : Ei) — > G{M ■ E2) commutes with the geometrical embeddings F E {Di) 
and F e (D 2 ). 

Thus for example if r\ + s\ = r 2 + s 2 then the Riemannian metric provides an 
isomorphism g : T^(M) —> T^(M) that pushes A£J to A^. Hence the corre- 
sponding functional calculus embedding commutes with the lowering or raising of 
indices. 
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